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Thanks to its intrinsic ability to preserve spin coherence, graphene is a prime material for spin¬ 
tronics. In this review article, we summarize recent achievements related to spintronics in graphene 
quantum dots and motivate this field from a spintronics and a materials science point of view. We 
focus on theory but also discuss recent experiments. The main sources of spin decoherence are inter¬ 
actions with lattice excitations and the hyperfme interaction with present nuclear spins. We explain 
effective spin-phonon coupling in detail and present a generic power law for the spin relaxation 
time T\ as a function of the magnetic field. For specific cases, we discuss spin relaxation in detail. 
The Heisenberg exchange interaction is paramount for coherent spin qubit operation and addressed 
in the context of magnetism in graphene nanoflakes. Nuclear spins in the host and surrounding 
material can be considered by several means and the influence of 13 C nuclei has been studied in 
detail. Impressive advances in general spintronics and the fabrication of graphene devices are likely 
to spark significant advances in spintronics with graphene quantum dots in the near future. 


I. INTRODUCTION 

Spintronics and graphene are exciting, vigorous, and 
rather new areas of research with a rapid pace of new 
discoveries. Books and a multitude of review articles are 
available for those that want to learn about one of these 
fields [1-9 . The intersection of both provides fertile soil 
for an abundance of fascinating physics and its fruits have 
been partially reaped by previous reviews [10 -12 . Here, 
we put the cherry on the cake and review specifically 
spintronics in graphene quantum dots. 

While classical electronics relies on the charge for infor¬ 
mation processing, spin-based electronics or spintronics 
is the paradigm of an advanced technology where the spin 
degree of freedom complements or even replaces charge 
as the carrier of information. The expected benefits of 
spintronics encompass non-volatile data storage, faster 
and more energy efficient data processing, increased data 
density, and many more. In order to fully exploit the po¬ 
tential, efficient generation, transport, transfer, manip¬ 
ulation, and detection of spin polarization is required. 
All these requirements are closely connected to material 
properties. Therefore, progress in spintronics seems in¬ 
trinsically connected to the quest for new materials with 
appropriate characteristics [1, 9, 12 . 

Graphene is a novel material that has undeniably 
stirred up the solid state community since its isolation 
more than a decade ago [13-16 . It is a monatomi- 
cally thin, quasi two-dimensional layer of carbon atoms 
arranged in a honeycomb lattice with two sublattices, 
A and B. Due to its atomic structure, it has the low¬ 
est surface mass density of all solids [17 . At the same 
time, its in-plane mechanical strength is second to no 
other material [18 . The hexagonal structure is the con¬ 
sequence of sp 2 -hybridized carbon atoms. The remain¬ 
ing p-orbitals form the 7r-bands responsible for transport. 
Electrons move quasi-relativistically with the Fermi ve¬ 
locity t;f~c/ 300, where c is the vacuum speed of light, 
and with high mobilities [19-23 . Spin-orbit couplings 
are small and nuclear spins are sparse, making it a po¬ 



FIG. 0: An electron spin in a graphene quantum dot interacts 
with the electron momentum via the spin-orbit coupling and 
with the sparse 13 C nuclear spin via the hyperfine interaction. 


tentially good host material for spin [24-29 . In addition 
to the real electron spin, the occupation of sublattices 
A and B can be described by the sublattice pseudospin. 
Moreover, in momentum space, the occupation of the 
two inequivalent Dirac points — where valence and con¬ 
duction bands touch — may be denoted by the valley 
isospin [4, 11 . Electron spin currents can be injected 
optically in graphene with Rashba-type spin-orbit inter¬ 
action [30-32J. This allows to leave away ferromagnetic 
contacts, thus eliminating a possible source of spin scat¬ 
tering. Magnetic behavior has been observed at the edges 
of graphene [33-35 . These and further properties make 
graphene a very promising material for spintronics. 

Quantum dots (QDs), in graphene as well as in other 
materials, are quasi zero-dimensional regions to which 
charge carriers can be confined [6 , 10 , 36 , 37 . With 
electrostatic gate electrodes that can be used to adjust 
the confinement potential and with electric contacts for 
transport, QDs allow in principle for full control over the 
individual electron. In particular, QDs provide a control¬ 
lable playground to investigate the behavior of spins as 
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well as sources of decoherence and methods to increase 
coherence times [7, 38 , 39 . In general, there are nuclear 
spins as well as electron spins in a solid. Unless explicitly 
stated otherwise, we always refer to electron spin. 

In this review, we focus on what we believe to be the 
most immediate questions related to spintronics with lo¬ 
calized electron spins in graphene. However, we point out 
that there are a several related systems that make use 
of carbon nanotubes or transition metal dichalcogenides 
(TMDCs) and where the electron spin degree of freedom 
might be substituted or complemented by a pseudospin 
or valley isospin [11, 40 -46 . The review is organized 
as follows. In the next section, we give an overview 
over spintronics in quantum dots and motivate the use 
of graphene. In Sec. Ill, we discuss various graphene 
quantum dots and their associated spin relaxation times. 
Closely related to the concept of QDs are localized states 
due to defects in graphene nanoflakes, which we treat 
in Sec. IV, In Sec. V, we address the influence of nu¬ 
clear spins on spintronics in graphene. Recent experi¬ 
mental achievements and the current status are the topic 
of Sec. VI, Finally, we give a summary and perspective 
in Sec. VII, 


II. SPINTRONICS 

Quantum effects become increasingly relevant as com¬ 
ponents of information processing devices shrink to the 
few-nanometer scale. This poses a challenge but also of¬ 
fers the chance to exploit quantum mechanics for data 
processing [47-49 . The spin is the canonical example of 
a quantum mechanical degree of freedom with two eigen¬ 
states. In solids, it is possible to use spin-orbit coupling 
and many other spin dependent effects to manipulate the 
spin with electric fields and hence much faster than with 
external magnetic fields. In order to prepare, manipu¬ 
late, and detect spin in an active and controlled way, 
spintronic devices may take advantage of any such effects 
[50-53. 

In contrast to electronics where carrier mobilities and 
lifetimes matter, spin mobilities and coherence times are 
relevant parameters for spintronics. The coherence of 
electron spins is mainly determined by the spin-orbit in¬ 
teraction (SOI), the abundance of nuclear spins in the 
host material, and magnetic behavior. These character¬ 
istics depend on the material and interlink spintronics 
with material science. Due to their lower magnetic mo¬ 
ment, nuclear spins couple weakly to the environment 
and exhibit typical coherence times orders of magnitude 
longer than for electron spins [1, 54 . 

Spintronics promises new applications and substan¬ 
tially better performance than standard micro- and na¬ 
noelectronics by exploiting the spin degree of freedom. 
Hard drive read heads are a multibillion dollar industry 
and magnetoresistive random-access memory (MRAM) 
is still waiting in the wings. Both technologies are based 
on the giant magnetoresistant effect, whose discovery in 


the late 1980s may be considered the birth of spintron¬ 
ics [55 -57 . Beyond classical computing, spintronics has 
the potential for key components of quantum computing 
[58, 59 . Quantum computing can be viewed as a coher¬ 
ent superposition of classical computations running in 
parallel, with algorithms like the one by Shor exploiting 
this [60 . Obviously, such algorithms require appropriate 
quantum hardware [61 . 

A. Spins in quantum dots 

As proposed by Loss and DiVincenzo, spin states of 
coupled QDs can be used to implement a universal set of 
one- and two-qubit gates (qubit is short for quantum bit). 
Spin qubit gates can be realized with magnetic fields and, 
using spin-dependent effects, also with electric fields, or 
a combination of both [58, 61. Such a solid state ap¬ 
proach to quantum computing benefits from the scal¬ 
ability known from standard electronics and from long 
coherence times [7, 62 -64 . 

The dot regions to which an integer number of elec¬ 
trons can be confined are typically on the length scale 
of 10-100 nanometers. The eigenstate spectrum is dis¬ 
crete and quantum dots are often referred to as artificial 
atoms. Coulomb repulsion and the exchange coupling 
also play a role [36 -38 . 

Also for spins that are confined to the QD, there re¬ 
mains a coupling to the environment. At finite tempera¬ 
tures, spins and phonons in the host material couple via 
the admixture mechanism, which involves spin-orbit in¬ 
teraction and electron-phonon coupling, or directly, via 
spin-orbit interaction only [65, 66 . Moreover, the hyper- 
fine interaction mediates a coupling between the electron 
spin and nuclear spins in its vicinity [38 , 67 . These mech¬ 
anisms lead to a finite spin lifetime, an important figure 
of merit for applications like the Loss-DiVincenzo quan¬ 
tum computer where it should exceed the clock time by 
orders of magnitude [61, 68 . The theoretical prediction 
and explanation of spin lifetimes (in graphene QDs) is 
part of this review. 


B. Why use graphene for spintronics? 

Its many outstanding properties make graphene, a 
monatomically thin layer of carbon atoms, an interesting 
material for a variety of applications, including spintron¬ 
ics. The low spin-orbit coupling and a small abundance of 
nuclear spins are beneficial for long spin coherence times 
[24 . In Table I, we compare the main sources of deco¬ 
herence for different materials. 

Electron mobilities in graphene are typically high and 
allow for spin transport over micrometer length scales 
[23 , 69 -72 . While graphene is naturally non-magnetic, it 
can exhibit magnetic behavior under certain conditions 
[33-35, 73 . More recently, there has been an effort to 
stack graphene and other two-dimensional materials like 
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hexagonal boron nitride (hBN) and TMDCs on top of 
each other, held together by van der Waals forces [74- 
77 . It is envisioned that by combining different materials 
in such van der Waals heterostructures, one may tailor 
physical properties. The concept is akin to bandgap en¬ 
gineering in semiconductor alloys and might once again 
highlight the connection between spintronics and mate¬ 
rials science. 

The hexagonal structure of graphene comes from the 
sp 2 -hybrid formed by three of the four outer-shell elec¬ 
trons in carbon. The resulting cr-bands are energetically 
far away from the charge neutrality point and play no role 
in transport. The remaining p z orbitals form the 7r-bands 
responsible for transport. In the low-energy approxima¬ 
tion, they obey the quasi-relativistic Dirac Hamiltonian 

[4, 

Ut, = hv F (ra x q x + a y q y ) + Aa z + V(x, y ), (1) 

where r labels the Dirac valleys K (r=+1) and K' 
(t=— 1) in reciprocal space, <j are Pauli matrices in stan¬ 
dard representation that describe the sublattice, and q 
is the wave vector measured w.r.t. K or K f , respectively. 
The sublattice splitting A vanishes for ideal graphene 
but may occur due to a substrate [76J. An electrostatic 
potential may be considered via the term V{x,y). The 
eigenspectrum of Eq. (1) with V =0 is given by 

E = ±yJ{hv F ) 2 (ql + ql) +A 2 (2) 

and has no gap if q =0 is possible and A=0. In that case, 
it is not possible to confine electrons electrostatically be¬ 
cause of Klein tunneling [78, 79 . In Sec. Ill, we will 
discuss special cases where the energy spectrum becomes 
gapped, thus enabling electrostatic confinement. 

Due to its low nuclear charge, carbon has a relatively 
low atomic spin-orbit interaction and this carries through 
to the spin-orbit interactions for band electrons in all car¬ 
bon based materials. Yet for flat graphene, the coupling 
is particularly weak because the i r- and cr-bands are or¬ 
thogonal [25 -28 . The spin-orbit coupling in graphene is 
given by 

Hsoi = Ai tg z s z + Ar (ra x s y - a y s x ), (3) 

where the first term denotes the intrinsic (or Dresselhaus- 
type) effect Hi and the second term stands for the extrin¬ 
sic (or Rashba-type) effect Hr, which may be induced 
by a substrate or an external electric field. The coupling 
strengths are not completely settled, but for concrete¬ 
ness we list Ai=12//eV and Ar=5 fieV x E[V/nm] (see 
Ref. [27 and also Table I). 

The two-dimensional surface states of three- 
dimensional topological insulators (TIs) are similar 
to graphene in that their low-energy excitations are also 
described by a Dirac Hamiltonian. While spin-orbit 
coupling is flimsy in graphene, it is strong in TIs. Future 
spintronics devices may thus rely on graphene for spin 
transport, combined with TIs for spin manipulation [11 . 


TABLE I: Sources of decoherence in different materials. If 
two spin-orbit splittings are listed, the upper (lower) number 
corresponds to the conduction (valence) band. Data on the 
different isotopes is taken from Ref. [29 . 
f Data for a (4,4) armchair carbon nanotube. The spin-orbit 
coupling scales inversely with the nanotube radius [80 . 
f The isotopes 92 Mo, 94 Mo, 96 Mo, 98 Mo, and 100 Mo have 
nuclear spin 0 and a combined natural abundance of 75%. 


Material 

Spin-orbit splitting Nuclear spin 

Graphene 

24 yeV [27 81. 

0 (99% 12 C) 

| (1% 13 C) 

Bilayer graphene 

24 /ieV [82 


Carbon nanotube^ 1.6 meV [12 


Silicene 

1.6 meV [83 

0 (95% 28,30 Si) 

| (5% 29 Si) 

Silicon (3D bulk) 

44 meV [84 


2D hexagonal BN 

15 /xeV 12 

| (80% n B) 


30 /ieV [12 

3 (20% 10 B) 


1 (>99% 14 N) 
I (<1% 15 N) 


2D MoS 2 

3meV 85 

0 (75% f Mo) 


147 meV 85 

| (25% 95,97 Mo) 



0 (99% 32 ’ 34 ’ 36 S) 



| (1% 33 S) 

GaAs (3D bulk) 

340 meV 84 

| (100% 69,71 Ga) 
f (100% 75 As) 


However, spin relaxation times measured in transport 
experiments in graphene are significantly shorter than 
expected and the underlying mechanisms remain some¬ 
what elusive [12, 69, 86J. Nevertheless, the past few 
years have brought many new insights and experimental 
spin lifetimes have increased with improvements in 
device fabrication [72, 87-90 . 

Since protons are 1836 times more massive than elec¬ 
trons, their magnetic moments are smaller than the Bohr 
magneton by that factor. This inhibits the exchange of 
angular momentum between electron spins and nuclear 
spins if a magnetic field is involved as the resulting Zee- 
man splitting scales linear with the magnetic moment. 
Still, nuclear spins do play an important role for spin co¬ 
herence [67 . The natural abundance of carbon isotopes 
is dominated by 99% 12 C, which has no nuclear spin. 
About 1% of carbon atoms are 13 C isotopes and have 
nuclear spin With the atomic distance of 1.42 A in 
graphene, one thus expects one atom with nuclear spin 
\ in a square area of (16 A) 2 . While such small QDs can 
be realized by electroburning [91^, typical QD dimensions 
are rather (16 nm) 2 or (160nm) 2 , according to 100 to 
10 000 atoms with nuclear spin [92 -96 . For silicon, whose 
natural abundance of isotopes is dominated by isotopes 
without nuclear spin (92% 28 Si, 3% 30 Si), electron spin 
coherence times exceeding seconds have been achieved af¬ 
ter isotopic purification [62 . Given that natural carbon 
has even fewer nuclear spins than natural silicon (see Ta- 
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ble I), this technique would be most suitable for graphene 
[97. 

In addition to the real spin, electrons in graphene have 
other binary degrees of freedom, namely the sublattice 
spin that refers to sublattices A and B (usually referred 
to as pseudospin) and the valley spin that refers to the 
Dirac points K and K' (technically also a pseudospm 
but usually called isospin). Moreover, the two layers in 
bilayer graphene can be attributed a layer spin. The 
valley spin also occurs in TMDCs and valleytronics - 
the electric initialization, manipulation, and detection of 
valley spin — is an active area of research [11, 44 , 45 , 98 . 
Now, we turn to graphene quantum dots and the lifetimes 
of electron spins confined to those QDs. 


III. GRAPHENE QUANTUM DOTS AND SPIN 
RELAXATION 

There are several possibilities to localize electrons to 
a well-defined area within a graphene sheet. One pos¬ 
sibility is to lithographically cut the desired shape of a 
graphene island into the 2D bulk, to which it remains 
connected by tiny nanoconstrictions. In this case, bound¬ 
ary effects need to be under control to understand and 
manipulate the behavior of electrons such QDs [99 , 100 . 
Other approaches involve inhomogeneous magnetic fields 
or rely on localization due to disorder [101, 102 . 

In order to benefit from the expertise with semicon¬ 
ductor QDs and to facilitate fast switching times, elec¬ 
trostatic confinement of electrons in gate-tunable QDs is 
desirable. For electron spins, the Dirac Hamiltonian im¬ 
plies two challenges, depicted in Fig. 1: first (i), the spec¬ 
trum, given by Eq. (2), must be gapped in order to avoid 
Klein tunneling, i.e., transmission of conduction electrons 
through the electrostatic barriers via the valence band 
states [79 . (ii) Due to spin and valley degeneracies, the 
orbital states of graphene electrons are fourfold degen¬ 
erate, i.e., they form four-level quantum systems rather 
than two-level qubits. Therefore, the time-reversal sym¬ 
metry between Dirac points K and K' needs to be lifted 
unless one accepts a complicated exchange interaction for 
the four-level system, which is also possible [98 . Both 
challenges are met by graphene nanoribbons with a cer¬ 
tain type of boundary condition [10, 24 as well as by 
circular quantum dots in mono- or bilayer graphene and 
with finite magnetic field [103 . 

In the following subsection, we introduce general con¬ 
cepts of spin relaxation via the spin-orbit interaction me¬ 
diated coupling to lattice vibrations. After that, we dis¬ 
cuss details of the mentioned graphene QDs with lifted 
valley degeneracies. 


A. Spin-orbit mediated spin relaxation 

If spin decoherence is dominated by interactions with 
phonons, nuclear spins can be neglected and the system 



FIG. 1: The electronic spectrum poses two challenges for spin- 
tronics with graphene quantum dots, (i) Pristine graphene 
has a linear, gapless dispersion (thin lines) that allows charge 
carriers to override electrostatic barriers. A finite energy gap 
is necessary to avoid this phenomenon known as the Klein 
paradox 78, 79 . (ii) The orbital states at valleys K and K' 
are connected by time reversal symmetry and hence twofold 
degenerate. Together with spin degeneracy, this valley de¬ 
generacy makes for a fourfold degenerate system rather than 
a two-level qubit, with an accordingly more complicated ex¬ 
change mechanism 98 . This can be avoided by lifting the 
valley degeneracy. If K and K' are degenerate due to time 
reversal symmetry, there will be two sets of Kramers qubits 
(solid and empty arrows, respectively), i.e., qubits that are 
connected by time reversal symmetry. Kramers spin qubits 
have long lifetimes due to the Van Vleck cancellation 65 , 104 . 


can be described by the Hamiltonian 


1~L — 77 e iec T 77phon T 77 S i 


K 


epc 5 


( 4 ) 


where H e iec and H p hon model the unperturbed electronic 
and vibrational systems, respectively [105 . The spin- 
orbit interaction 'H SO i= / H Y +'H[ is the sum of extrinsic 
(or Rashba-type) coupling and intrinsic (or Dresselhaus- 
type) coupling. Together with the electron-phonon cou¬ 
pling Hepc, h mediates the coupling between spins and 
the lattice phonons. We denote the orbital eigenstates 
of Heiec as | n). For graphene, l~i e i ec is usually given by 
some form of Eq. (1). The pure vibrational modes are 
described by 


^phon — ^ ^ 


a,q 




( 5 ) 


where the summation runs over all acoustic phonon 
branches a and wave vectors q. The angular frequency 
uj a , q of a harmonic vibrational mode is implicitly deter¬ 
mined by a and qr, and m a ^ q is the occupation number 
of this mode. The according eigenstates are determined 
by these occupation numbers, \m a , q ). 

Phonons are bosons with zero spin and hence cannot 
couple to electron spin directly but via the admixture 
mechanism which combines Rashba-type spin-orbit inter¬ 
action H r and electron-phonon coupling % epc [65 . Let 
|n, s}(°) be the product state of an electron in the orbital 
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TABLE II: Overview of the parameters used for the generic 
model of T^ 1 (B) in Eq. (9). 


Description 

Parameter Value 

Dipole approximation 

a 

0 (zeroth order) 

1 (first order) 

Kramers pair 

p 

0 (no) 

1 (yes) 

System dimension 

D 

1 

o 



z 

3 

Phonon dispersion 

g 

1 (linear) 

2 (quadratic) 

... (higher) 


state | n) and with spin |s). In lowest order, H r admixes 
product states with opposite spin, 


|nt) = |nt) (0) + £ l n 'E 0) 

n'^n 


(°V + |^ r |nt) (0) 

B n En'+g^BB 


(6) 


where E n is the orbital energy and EgpsB is the Zeeman 
energy with the ^-factor (g=2 for graphene [106 ), the 
Bohr magneton //p, and a magnetic field B. An according 
expression for | n leads to finite spin-flipping matrix 
elements (n i \H eY)C \n t) = 


E 

n'^n 


(H epc ) nn '(H I )$ n | (^epc)n>n(K)t;T 
B n E n r+gfiBB E n E n r ggsB J 


( 7 ) 


where we denote the numerator in Eq. (6) as (% r )^t n and 
the spin-conserving transitions of (H epc )n , n accordingly. 
The working principle of the admixture mechanism is il¬ 
lustrated in Fig. 2 (a). With p(hu: a ^ q ) as the phonon 
density of states, the relaxation rate can now be calcu¬ 
lated via Fermi’s golden rule, 

2^ _ 

= _ ^T" ^ > I ^ ^o:,q~l~l \B e p C \fl , TTi aq ) | p{7l^a,q) • 
a,q 

( 8 ) 

Due to energy conservation, the phonon energy must 
match the Zeeman splitting, Tiuj a , q =gpBB. For typ¬ 
ical laboratory magnetic fields, E><30T, this implies 
long-wavelength acoustic phonons at the center of the 
Brillouin zone that can be described by the continuum 
model. In a very general picture, the atomic displace¬ 
ment in the presence of an acoustic phonon behaves as 

_ i 

'Ua,q{r)(xuj a , q e' l ( q ' r - UJcx ’ qt \ where the prefactor occurs in 
the normal coordinate for any harmonic oscillator. A 
constant acoustic displacement u corresponds to a mere 
translation of the lattice and hence does not give rise to 
electron-phonon coupling. Usually, H epc involves deriva¬ 
tives diUj, see, e.g., Eqs. (10) and (11), or Ref. [65 . 

The matrix element given by Eq. (7) scales with 
q 1+a B^~^, w here one power in q comes from these spa¬ 
tial derivatives. In the long wavelength limit, the fac¬ 
tor e iqr can be replaced with the dipole approximation 


TABLE III: Specific examples for Eq. (9). 
f In (armchair) graphene nanoribbons, 'Hdef and Hblc cor¬ 
respond to the same parameters and Hdsp does not give a 


contribution, see Subsec. 

III B 

or 

Ref. 

[109 . 

System 

a 

p 

D 

v B 

^ocB- 

Graphene nanoribbon^ 

1 

1 

1 

1 5 

[109 

Bulk graphene ('Hdef) 

1 

0 

2 

1 4 

[110. 

Bulk graphene (Hblc) 

0 

0 

2 

1 2 

[110. 

Bulk graphene (Hdsp) 

0 

0 

2 

2 0 

[110. 

Bulk graphene (Hdsp) 

0 

0 

2 

1 2 

[110. 

GaAs 

0 

1 

3 

1 5 

[65, 111 

Silicon metal-oxide 

1 

1 

3 

1 7 

[112. 


e iq r ~l+iq r. Depending on the symmetry of the elec¬ 
tronic states, the constant term may (<a= 0 ) or may not 
(a= 1) give a contribution to the matrix element. If the 
two states of the spin qubit are related by time rever¬ 
sal symmetry in the zero field limit, i.e., if they form a 
Kramers pair (Fig. 1 ), the matrix element will vanish un¬ 
less the Zeeman energy in the denominators is taken into 
account [65 , 104 . Therefore, /3=1 (/3=0) if the spin states 
form a (no) Kramers pair. Finally, the factor B~ 2 comes 
from the normal coordinate of the acoustic phonon. 

The density of states in Eq. ( 8 ) scales with where 
Aocg D_1 is the content of the iso-energy hyperplane in D- 
dimensional reciprocal space (De{ 1,2,3}). Assuming a 
phonon dispersion ujcc q 71 ( 77 EN), the derivative scales as 

1 — V D 1 

uj f and hence pocB n , where we have used uxxB. Us¬ 
ing qocBv for the matrix element, too, we find a general 
estimate for Eq. ( 8 ), 

7 , 1 x) Bn . (9) 

An overview of the parameters in this formula is pro¬ 
vided in Table II and specific examples are shown in 
Table III, Unless direct restoring forces are absent (as, 
e.g., for out-of-plane modes in free graphene [107, 108 ) 
acoustic modes disperse linearly ( 77 = 1 ) and the above re¬ 
lation simplifies to T-f X (xi? 2 ( a+ ^ +jD . Eq. (9) suggests a 
monotonic behavior of T 1 _ 1 (,B). For large enough mag¬ 
netic fields, non-monotonic behavior may occur due to 
(possibly avoided) crossings of orbital levels, destructive 
interference of different mechanisms, or Van Hove singu¬ 
larities in the phonon density of states, see e.g. Fig. 4 
[109 , 113 -116 . 

In graphene, we consider two mechanisms for the 
electron-phonon coupling, namely the deformation po¬ 
tential 'HDEF=gi(u xx +Uyy) and the bond-length change 

' H B LC=T(T x (TA x )+(Ty(TAy ), 

%EPC='%DEF+'%BLC • ( 10 ) 

In a minimal coupling picture, the form of Hblc, where 
(A x ,Ay)=g 2 (u xx -Uyy,-2u X y), highlights the equiva- 
lence of strain and a valley-dependent magnetic field, 
see Eq. (1) and Refs. [117, 118 . The strain tensor is 
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(b) Z 

s n 




FIG. 2: (a) The admixture mechanism allows an electron in 
orbital state | n) to flip its spin via a virtual transition to a 
different orbital state | n). Both electron-phonon coupling 
Tie P c and Rashba-type spin-orbit interaction Ti r couple the 
different orbits but only Ti r leads to a spin flip, thus effectively 
enabling a spin flip within the orbital state | n). (b) Due 

to out-of-plane phonons, the local coordinate frame of the 
graphene lattice, £', that describes pseudospin cr' may differ 
from the laboratory frame £ that fixes the quantization axis 
of the electron spin s in an external magnetic held. In the 
local lattice frame, s z may depend on all components of s, as 
indicated by the red arrows. This changes the intrinsic spin- 
orbit coupling to the form of Eq. (11), where the primes have 
been dropped. 


given by Uij=(diUj+djUi)/2 and the coupling constants 
are gi«30eV and #2~l-5 eV [110, 117 . We point out 
that all possible in-plane strains Uij (i,j=1 ,2) are con¬ 
sidered by Eq. (10). 

The admixture mechanism creates a lowest-order con¬ 
tribution for in-plane deformations ( u x ,u y ) but not for 
out-of-plane deformations u z since a local tilt of the 
graphene lattice does not give rise to finite Tippc • Out-of- 
plane phonons can still induce a spin flip via the intrinsic 
spin-orbit interaction Ti \, given in Eq. (3), or in higher 
order, where stretching effects lead to nonzero Hepc- 
In the former case a small local tilt ( d x ,d y )u z , where 
\{d x ,d y )u z |<Cl, changes the mutual orientation of pseu¬ 
dospin <j and real spin s. Asa consequence, the intrinsic 
spin-orbit coupling acquires spin-flip terms, 

'Hi=^zTz(s z -d x u z s x -dyU z s y )='Hi + 'Hdsp , (11) 

and thus gives rise to direct spin-phonon coupling Tipsp 
[66 . The origin of this mechanism is illustrated in Fig. 2 
(b). Below, we discuss the relaxation times of the elec¬ 
tron spin in specific graphene QDs where the electron 
spectrum is gapped and where the valley degeneracy is 
lifted. 


B. Armchair graphene nanoribbons 

In graphene nanoribbons (GNRs), electrons are natu¬ 
rally confined within the quasi one-dimensional structure 
and can move freely only along the longitudinal direction 
(y -axis in Fig. 3). 

While the longitudinal wavenumber q y may vary con¬ 
tinuously, the transverse wavenumber q n is discrete and 
depends on the ribbon width W as well as the specific 
boundary conditions. In GNRs with lateral armchair 



FIG. 3: Sketch of (a) an armchair GNR and (b) its electronic 
spectrum, (a) Here, the width of the GNR is characterized by 
two integers, m— 3 and /jl= — 1. The resulting bandgap allows 
for electrostatic confinement along the y- axis. A potential 
V(y) defines the barrier regions B1 and B2 (shaded) as well 
as the dot region D in between, (b) For A V<E gap , there is 
one or more bound state(s) inside the quantum dot (green 
lines). A continuum of extended states exists for energies 
above the confinement potential (pink shade). Both bound 
and extended states provide virtual intermediate levels for 
the admixture mechanism responsible for spin relaxation. 


boundaries, all components of the wavefunction spinor 
must vanish on the outermost atoms, which implies 

q n =ir(n-n/3)/W, (12) 

where 77=0, ±1, ±2,... labels the transverse excitations. 
The width of the armchair GNR, W=(3m-\-/jL)y/3a, is de¬ 
termined by ra=l,2,3,... and /iE{ —1, 0,+1}, see Fig. 3 
(a) and Ref. [119 . For a ribbon with fi=± 1, all transverse 
wavenumbers q n will be finite, thus leading to a gapped 
phonon spectrum with E gav =2TiVY\qo\=2i:TiVY/3W . For 
W=30nm, this amounts to £ g ap~40meV. The gap al¬ 
lows for electrostatic confinement of electrons along the 
longitudinal direction by virtue of an electrostatic po¬ 
tential V(y). Analytical solutions for bound states are 
known if the potential has the piecewise constant form 
V(y)=AV(0(—y) + 0(y—L )), where L is the longitudi¬ 
nal dimension of the quantum dot, see Fig. 3 (a) and 
Ref. [24 . At least one bound state exists if the confine¬ 
ment potential does not exceed the bandgap, A V<E gap . 
Depending on AV, there may be several bound states 
and independent of AV, there is always a continuum of 
extended states energetically above AV. Since £ gap is 
typically of the order of AV, we consider both bound 
and extended states for virtual transitions within the ad¬ 
mixture mechanism. 

Due to the boundary conditions, the electronic states 
become superpositions of K and K' with equal weight 
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FIG. 4: Spin relaxation rate T^ 1 (B) of the QD ground state 
for L/W— 5, AV=1.8hvFqo, and free mechanical boundaries 
[109 . The doubly logarithmic scale highlights the behavior 
T^ 1 cxB 5 for cZKO. 5, in agreement with Eq. (9). The dimen¬ 
sionless phonon frequency O is explained in the main text. 
Destructive interference of processes originating from 7-Ldef 
and LLbec causes two dips of T x _1 for 2<0<3 and a Van Hove 
singularity of the phonon density of states leads to a diver¬ 
gence around 0= 3. 


on both valleys. This has important consequences. First, 
the valley degeneracy is lifted, thus circumventing the ne¬ 
cessity of a complicated exchange interaction. Neverthe¬ 
less, the orbital system is symmetric under time reversal 
and the spin states in a certain orbit form a Kramers pair 
at zero magnetic field, which suppresses spin relaxation 
due to Van Vleck cancellation. Furthermore, the direct 
spin-phonon coupling %dsp does not allow for spin relax¬ 
ation as the according matrix element is proportional to 
the valley index r and hence vanishes after summing up 
the contributions from K and K'. The spin may still re¬ 
lax via Repc and the admixture mechanism. In Eq. (8), 
phonon emission [absorption] is proportional to m a , q +1 
where the phonon number m a , q is determined by 
the Bose-Einstein distribution. Here, we assume phonon 
vacuum, i.e., pure phonon emission. 

The spin relaxation depends on the confinement po¬ 
tential AF, the QD aspect ration L/W , and the applied 
magnetic field T>, which is assumed to be perpendicular 
to the graphene plane. Moreover, the mechanical bound¬ 
ary conditions determine the phonon spectrum and hence 
also the relaxation rate. In Fig. 4, we plot the spin re¬ 
laxation rate T^/ 1 (B) in blue as calculated with Eq. (8) 
for L/W= 5, AV=1.8hvEqo : and free lateral boundaries 
[108, 109 . In agreement with Eq. (9), the doubly loga¬ 
rithmic scale highlights the behavior T 1 _1 oc5 5 for u)<0.5, 
where 0=uj^/ p g /YW is the dimensionless frequency of 
in-plane phonons with Y and p g as Young’s modulus 
and the mass density of graphene, respectively [108 , 109 . 
Due to destructive interference of processes originating 
from 'Hbef and Relc, the relaxation rate features two 
dips for 2<c9<3. A Van Hove singularity of the phonon 
density of states p leads to a divergence around 0=3 — 
a common effect for quasi one-dimensional systems [116 . 

Results for different values of AV and L/W , or for 
fixed mechanical boundary conditions can be found in 


Ref. [109 . As one would expect, virtual transitions to ex¬ 
tended states gain importance as AV/E gap decreases and 
even dominate over virtual transitions to bound states 
for AV/E’gap—^O. Fixed mechanical boundaries lead to a 
gapped phonon spectrum which excludes spin relaxation 
via single-phonon processes for O<2.06 (or T><8.25T for 
IF=30 nm). 

As a quasi one-dimensional system, the armchair GNR 
is particularly interesting for spintronics-based quantum 
computing as it leaves two spatial dimensions for sys¬ 
tem periphery in a scaled-up device with multiple QDs 
along the GNR. The wave functions of bound states also 
extend to neighboring QDs and thus allow for multiple- 
qubit operations. Moreover, the small bandgap allows 
the coupling of non-neighboring QDs via detuning of the 
intermediate QDs, which allows for more efficient quan¬ 
tum computation since shuttling of qubits can be avoided 
and thus raises the threshold for fault-tolerant quantum 
computing [10, 24, 120, 121 . 


C. Mono- and bilayer graphene disks 

Gate-tunable QDs with an electrostatic confinement 
potential can exist in 2D bulk mono- and bilayer 
graphene. Analytical solutions exist if the systems have 
circular symmetry in the graphene plane and the confine¬ 
ment potential has a rectangular profile, V(r)=AV6(r — 
R), where r is the radial coordinate and R is the radius 
of the disk-shaped quantum dot. Due to the absence 
of intervalley scattering, the orbital degeneracy between 
valleys K and K' can be broken in a controlled way by 
the inclusion of a mass term — a sublattice splitting for 
monolayer graphene or a layer bias in the case of bilayer 
graphene — in combination with a magnetic field. In 
particular, the splitting of orbits in opposite valleys can 
be tuned with the applied magnetic field. In addition to 
a valley splitting, the mass term also gives rise to the 
required bandgap [10, 103]. For magnetic fields, where 
the magnetic length lE=\/Ti/eB is much smaller than R, 
the discrete QD spectrum converges to the bulk Landau 
levels. Below, we outline the analytical solution of the 
eigensystem both for mono- and for bilayer graphene. We 
first treat monolayer graphene disks and bilayer graphene 
afterwards. 

As is obvious from Eq. (2), a sublattice potential A 
will open up a bandgap of E gap =2 A. A sublattice po¬ 
tential arises when graphene lies on a substrate but the 
two sublattices couple differently to the substrate ma¬ 
terial as, e.g., for graphene on hBN, where A>100meV 
has been observed [76 . While electron spin is neglected, 
here, the perpendicular component of the magnetic field 
is taken into account via minimal coupling to the motion 
of the charge carriers. Pseudospin can be included in the 
total angular momentum, J z =-id^-\-a z /2^ which com¬ 
mutes with the Hamiltonian. As a result of the circular 
symmetry, the azimuthal dependence can be separated 
away, leaving a differential equation for the radial depen- 
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FIG. 5: Spin relaxation times Ti(5) for the ground state 
of a circular QD in monolayer graphene and a perpendicu¬ 
lar magnetic held. The QD radius is R =25 nm and here, 
the bandgap 5 gap =2 A is twice the confinement potential, 
AV=A=260 meV. The coupling constants gi and #2 indi¬ 
cate the associated mechanisms, 5def and 5blc, respec¬ 
tively. Both LA and TA in-plane modes disperse linearly 
while the out-of-plane mode has a transitions from quadratic 
(ZA) to linear (ZA’) dispersion. The sum of all relaxation 
channels is represented by the red dotted, blue dashed, and 
black solid lines. A quadratic (linear) dispersion of the out- 
of-plane modes is assumed for the red (blue) curve, and a 
crossover between them for the black curve. The inset shows 
the universal dependence of T^ 1 ocf(6) on the inclination an¬ 
gle 0 of the magnetic held [110 . ©APS. 


dence. Both for the dot region ( r<R ) and for the barrier 
region (r>5), the solutions are found to be hypergeo- 
metric functions. The matching condition for r=R leads 
to a discrete spectrum of bound states. The valley de¬ 
generacy of K and K' is lifted for hnite magnetic held 
and exceeds the Zeeman splitting for experimentally ac¬ 
cessible parameters [103 . 

As for spins in armchair GNR quantum dots, relax¬ 
ation via the admixture mechanism and via direct spin- 
phonon coupling will be considered in the following. As¬ 
suming A=260meV, no continuum of extended states 
has been considered in Ref. [110 . As a result of the 
broken valley degeneracy, the spin states do not form a 
Kramers pair and Van Vleck cancellation does not occur 
in circular QDs in monolayer graphene. The spectrum 
of acoustic phonons consists of three branches, trans¬ 
verse (TA) and longitudinal (LA) in-plane modes, and 
one out-of-plane mode (ZA). Near the center of the Bril- 
louin zone, the in-plane modes disperse linearly and the 
ZA modes quadratically with a crossover to linear be¬ 
havior. Only the LA mode changes the size of the unit 
cell and is thus the only mode where the deformation 
potential (5def) is nonzero. Both LA and TA deform 
the unit cell and hence generate a finite coupling via the 
bond-length change (5blc)- ZA modes give rise to direct 
spin-phonon coupling (5dsp)- The according spin relax¬ 
ation times for an electron in the QD ground state and 


phonon vacuum have been calculated in Ref. [110 with 
specific parameters. The results are shown in Fig. 5. For 
excited QD states, the relaxation times decrease quickly. 

In a two-dimensional material, the minimal coupling 
prescription involves the perpendicular components of 
the 5-field. But the Zeeman energy, which is trans¬ 
formed to lattice excitations by the mentioned couplings, 
depends on the entire magnetic field strength. As a re¬ 
sult, all relaxation rates also depend on the inclination 
angle 0 between the vector of the magnetic field and the 
normal vector of the graphene plane. For purely in-plane 
magnetic fields, the relaxation rates I© 1 are reduced by 
50% compared with the case of a perpendicular field, see 
inset in Fig. 5, 

In bilayer graphene, a gap opens if the electrostatic 
potential is different for the two layers. Then, an electro¬ 
static confinement potential can be used to confine charge 
carriers in a quantum dot. Such systems are studied in 
experiment [93-96 . With appropriate top and bottom 
gates, the bandgap as well as the confinement potential 
can be tuned electrically and independently [103 . In a 
simple description, the system consists of two individual 
graphene layers, where one sublattice of the upper layer 
couples to one sublattice of the lower layer with an inter¬ 
layer hopping amplitude £i2~0.4eV [122 . The inversion 
symmetry of the bilayer system is broken by the potential 
difference between the distinct layers. 

As for the case of monolayer graphene above, the prob¬ 
lem can be solved analytically for a potential of the form 
V(r)=AV0(r — R ). Due to the contribution of the pseu¬ 
dospin of both layers, the total angular momentum is an 
integer here. The analytic solution of the eigensystem 
follows the same procedure as for the monolayer case but 
is more complicated because of the two coupled layers. 
The radial parts of the wavefunctions again turn out to 
be hypergeometric functions both inside and outside the 
dot region, and the matching condition at r=R leads to a 
discrete spectrum of bound states. The valley degeneracy 
for 5=0 can be lifted with a finite magnetic field perpen¬ 
dicular to the graphene plane and for Ib^R, the energy 
spectrum approaches the Landau level structure. In par¬ 
ticular, there are no Kramers pairs for a finite magnetic 
field perpendicular to the graphene plane. 

To our knowledge, the spin relaxation time T\ of a 
bound state inside the QD has not been estimated yet. 
For acoustic phonons neat the center of the Brillouin 
zone, the atoms within a unit cell do not vibrate against 
each other. Moreover, the interlayer Van der Waals 
bonds of bilayer graphene are much less rigid than the 
intralayer cr-bond. As a consequence, the change of the 
interlayer distance can be neglected and the individual 
layers can be treated completely independently for small 
wavenumbers. The deformation potential and the bond- 
length change within each layer are thus given by two 
copies of Eq. (10), one for each layer. Using the spin- 
orbit coupling of bilayer graphene, T-f 1 can in principle 
be calculated via the admixture mechanism as for mono- 
layer graphene [82 , 123 . As for monolayer graphene, 
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out-of-plane modes change the mutual orientation of real 
spin s and lattice related spins (pseudospin a and for bi¬ 
layer graphene also the layer spin ^ [82 ). The resulting 
direct spin-phonon coupling might turn out to be more 
complicated than that for monolayer graphene. For cou¬ 
pling to the in-plane modes of bilayer graphene, Eqs. (9) 
and (10) suggest a behavior T^ 1 ocB 2 or ocF> 4 , depending 
on the order of the dipole approximation for the specific 
mechanism. 


IV. MAGNETISM IN GRAPHENE 
NANOFLAKES 


Both ferromagnetic and antiferromagnetic materials 
are important for data processing and storage (e.g. 
MRAM) |1 . In contrast to typical magnetic elements 
like Fe, Ni, Co, and rare earths, carbon does not con¬ 
tain occupied d- or f-orbitals and pristine graphene is 
strongly diamagnetic [34J. Nevertheless, the prospects of 
magnetic ordering at room temperature and unmatched 
data density in monatomically thin layers motivate in¬ 
tense research efforts aimed at magnetism in graphene. 
Magnetic moments may appear in graphene because of 
vacancies [33 , light or heavy adatoms [124 -126 , edge 
effects [35 , 127 , or molecular doping [128, 129J. The 
resulting magnetic structures may be 0D, ID, or 2D, de¬ 
pending on the origin of magnetism (e.g., isolated de¬ 
fects, GNR edges, or molecular doping of bulk graphene, 
respectively). The experimental detection of magnetic 
moments is possible via SQUID (superconducting quan¬ 
tum interference device) magnetometry [124 , spin trans¬ 
port measurement [125J, or spin-sensitive STM (scanning 
tunneling microscopy) [128, 130 . 

In the following, we focus on magnetic moments in¬ 
duced by vacancies or light adatoms. According to Lieb’s 
theorem, the magnetic moment of the ground state of the 
Hubbard model for graphene is given by /i b\Na — 1Vb|, 
where JVa and Nb are the numbers of sublattice sites 
[131, 132 . Carbon atoms can be removed from the 
graphene lattice by means of irradiation with electrons or 
argon ions [125, 133, 134 . Alternatively, light adatoms 
like H and F passivate the p z orbital of the carbon atom 
to which they bond, thus effectively removing this car¬ 
bon atom from its sublattice. A hBN substrate may be 
used to stabilize the hydrogen adsorption on one sub¬ 
lattice and to suppress migration of the adatoms [135 . 
Curie temperatures for the ferromagnetic ground state 
can exceed 300 K [33 , 136 . 

The above defects, i.e., vacancies and light adatoms, 
are similar to quantum dots in that they lead to localized 
states. In Ref. [73 , we have studied the exchange cou¬ 
pling between two spin states localized at two separate 
yet nearby vacancies in a graphene nanoflake (GNF). The 
need for a bandgap does not play a role in GNFs since 
electrons are naturally confined within such quasi zero¬ 
dimensional structures. Specifically, we consider GNFs 
with hexagonal symmetry — as found for flakes grown 
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FIG. 6: On-site probability densities for an armchair termi¬ 
nated hexagonal graphene nanoflake with a total 682 sites and 
two vacancies at r vac =(0, ±8a). (a) Antibonding eigenstate 
|1). Here, we only consider nearest neighbor hopping such 
that the (anti-)bonding energy eigenstates 1341) and 1342) 
lie in the middle of the spectrum and have the same on-site 
probability densities, (b) localized vacancy state | -\-y). The 
probability density of |— y) looks similar yet mirrored about 
the x-axis 73 . 


by chemical vapor deposition (CVD) [137, 138 — and 
either zigzag or armchair terminations. Two atoms are 
removed at sites r vac =(0, ±7/) w.r.t. the flake center such 
that the lattice retains the mirror symmetries A4 X :x\-^—x 
and M y :y^—y, see Fig. 6. 

We describe the electronic system with a full tight- 
binding calculation that involves hoppings t^{B) be¬ 
tween lattice sites i and j up to third nearest neighbors 
(n= 1, 2, 3). A perpendicular magnetic field is included in 
the hoppings via the Peierls phase. The resulting spec¬ 
trum is discrete and may or may not feature degeneracies, 
depending on the exact GNF configuration, H-field, and 
maximal order of nearest neighbor hoppings. 

The two nearby vacancies can be interpreted as a dou¬ 
ble quantum dot (DQD), whose localized states form a 
bonding and an antibonding eigenstate within the GNF 
spectrum. If the energy splitting A between the bonding 
and antibonding eigenstates is much smaller than their 
energy differences to the remaining spectrum, the orbital 
part of the DQD can be described by 

^ DQD = f) , (13) 

where we employ the localized states {|-by), | —y)} as a 
basis, E is the mean energy of the (anti-)bonding states, 
and where the intervacancy hopping with \t\=A/2 may 
acquire a Peierls phase due to the magnetic field. Sym¬ 
metry allows us to find the explicit form of |=b y), see 
Fig. 6. We complement this picture with an on-site 
Coulomb repulsion U and aim to calculate the exchange 
coupling J, which energetically splits the singlet state 
|S') = ^(c 1 f2/t c^^-c 1 f ^c^ t )|0) from the triplet state 

l T °^ = ^^ c +yt c -yt +c +?/t c -yt)l 0 ^ and determines their 
mutual dynamics. Within this model and for |£|<CU, we 
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FIG. 7: The exchange coupling J (solid red line and left axis) 
and the eigenenergies (dashed black lines and right axis) of 
the according (anti-)bonding states are plotted against the 
perpendicular magnetic field. The shown data belongs to a 
nanoflake with zigzag terminations and a total 598 atoms 73 . 
Two vacancies are located at r vac =(0, ±10a). Only nearest 
neighbor hopping is taken into account, such that the (anti- 
bonding energy eigenstates 1299) and 1300) lie in the middle 
of the spectrum. One flux quantum passes through the flake 
at 5=145.8 T. 


find that the singlet state is preferred, i.e., antiferromag- 
netic ordering, and 


J= 


4 |*£ 

U 


(14) 


for the exchange coupling. Both \t\ and 7/, and hence 
J, depend on the perpendicular magnetic field that is 
included from the beginning. Typically, the magnetic 
field has a significant effect when the total magnetic flux 
through the GNF, 4>=A5, where A is the surface area 
of the nanoflake, reaches one flux quantum, <F 0 =/i/2e, 
where h is Planck’s constant and e is the unit charge. 

Depending on the overall size of the GNF, the distance 
between the vacancies, and the GNF termination (zigzag 
or armchair), the 5-field dependence can be used to tune 
J over several orders of magnitude. For some flake con¬ 
figurations, it is possible to tune the spectrum into a 
degeneracy, A(5)— t(B)=J(B)=0, thus switching off the 
antiferromagnetic order, see Fig. 7. Such an in-situ tun- 
ability of the exchange coupling is most promising for 
spintronics applications as it allows to change the mag¬ 
netic behavior without preparing a new device. For flakes 
with armchair edges, we find that J increases for GNFs 
of larger overall size and for smaller separation distances 
of the two vacancies. Bulk graphene can be viewed as 
the limit where the GNF diameter is much greater than 
the vacancy separation and where edge effects become 
negligible. We expect that the former trend saturates in 
this limit. The latter trend applies separately for two 
cases: (i) when the atom that has been removed for the 
vacancy at (0, +y) belongs the A sublattice and (ii) when 
it belongs to the B sublattice. Such trends are not obvi¬ 
ous for zigzag terminated flakes, where edge effects play 
a significant role for GNFs with up to ^10 4 atoms [73 . 
Neel temperatures, Tn=J/£)b with the Boltzmann con¬ 
stant &b 5 range from below 4 K to beyond 300 K and thus 
allow for experimental verification of these results, possi¬ 


bly via the techniques mentioned at the beginning of this 
section. Including a non-local Coulomb interaction might 
lead to ferromagnetic ordering, J<0, and possibly enable 
an in-situ tunable crossover from antiferromagnetism to 
ferromagnetism [139 . 

Our method relies on symmetry and hence the removal 
of an A and a B atom. Nevertheless, ferromagnetic or¬ 
dering is expected if the two atoms are removed from 
the same sublattice and only on-site Coulomb interac¬ 
tion is considered [34, 140 . We also point out that the 
exchange coupling J can be used for coherent two-qubit 
operations, as required for universal quantum computing 
[39, 58, 63 , 139 . 


V. NUCLEAR SPINS AND GRAPHENE 

Naturally, only 1% of all carbon atoms belong to the 
isotope 13 C, which is the only stable carbon isotope with 
finite nuclear spin, see Table I This abundance can be 
lowered by isotopic purification [97 . Possibly, this is 
the reason why the combination of nuclear spins and 
graphene has so far received little attention from the sci¬ 
entific community. However, there are at least three sit¬ 
uations where nuclear spins become important for spin¬ 
tronics in graphene, (i) First, it has been found the¬ 
oretically that T\ relaxation times of electron spins in 
graphene QDs can exceed 100 /rs for certain parameters 
and if nuclear spins are neglected, see Figs. 4 and 5 or 
Refs. [109, 110 . In this case, the spin dephasing time T\ 
might be limited by the coupling to nuclear spins of the 
13 C nuclei, (ii) Then, state-of-the-art graphene QDs of¬ 
ten employ a graphene (bi-)layer encapsulated by bottom 
and top sheets of hexagonal boron nitride. Neither boron 
nor nitrogen posses any stable isotopes without nuclear 
spin, see Table I. The spin of an electron in such an en¬ 
capsulated graphene QD is in immediate contact with the 
nuclear spins of hBN. (iii) Finally, the nuclear spin of an 
impurity, of an adatom, or of 13 C itself might be interest¬ 
ing for spintronics. Natural monolayer graphene contains 
one nuclear spin per (16 A ) 2 surface area. This density 
can be controlled by enriching or depleting 13 C. Possibly, 
control and coupling of nuclear spins in the graphene 
sheet or on top of it can be achieved. Ultimately, this 
could lead to 2D nuclear spintronics with similarities to 
nitrogen vacancy centers in diamond, where 13 C nuclear 
spins have been used for quantum operations [141, 142 . 

The Hamiltonian for an ensemble of nuclear spins in a 
solid has five different contributions [67 , 

l^nuc—5 Z + 5hf + 5orb + 'Hdd + 5 q . (15) 

The term 7~L Z arises due to the Zeeman energies of nuclear 
spins in a magnetic field B. The second term describes 
the hyperfine interaction between electron spins and nu¬ 
clear spins. It is the sum of two terms, 'Hhf='H c +'H a . The 
contact (isotropic) hyperfine interaction can be modeled 
as an Overhauser field, 5 n, that acts on the electron spin 
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s in a similar way as an external magnetic field, 


H c = A k I k J • s=g/i B #N • s . (16) 

Here, I k is the spin operator for the nucleus at site k 
and the coupling strength A k is proportional to the lo¬ 
cal probability density of the electron and the gyromag- 
netic ratio of the nucleus at site k [38, 67, 143, 144 . 
Eq. ( 16) vanishes for 7r-band electrons (in particular those 
of graphene) as their probability densities at the nuclear 
sites are zero due to odd symmetry. In this case (in 
particular for graphene), the anisotropic hyperfine inter¬ 
action, H a , and the coupling of nuclear spin and electron 
orbit, H 0 rb, dominate electron-nuclear coupling. In ad¬ 
dition to Hhf, also H 0 rb can mediate a coupling of nu¬ 
clear and electron spins if spin-orbit interaction is in¬ 
cluded. The nuclear dipole-dipole Hamiltonian, Had, di¬ 
rectly couples the magnetic moments of distinct nuclear 
spins. And nuclei with spin />| (thus excluding 13 C) 
interact with an electric field gradient via the nuclear 
quadrupolar coupling, H q . All these terms are discussed 
in detail in Ref. [67 . 

Spintronics with nuclear spins — point (iii) in our list 
above — would rely on the coupling of nuclear spins to 
external fields, which is achieved by H z and H q , as well 
as the coupling between different nuclear spins, which 
is mediated by Had- To our knowledge, the effect of nu¬ 
clear spins from hBN on the electron spin in encapsulated 
graphene QDs — point (ii) — have not been studied, 
yet. In the remainder of this section, we review studies 
on point (i), the decoherence of the electron spin due to 
hyperfine interaction with the nuclear spins of 13 C. 

The contact hyperfine interaction vanishes for flat 
graphene and indirect coupling of electron and nuclear 
spins via H or b and spin-orbit interaction is expected to 
be small and hence neglected. Consequently, the cou¬ 
pling between electron and nuclear spins in graphene is 
given by 

^hf- ^ ^ A k x I k x S X -\-A k yI k ySy-\-A k z I k z S Z , (17) 

k 

where k indexes the sites of 13 C isotopes. The coupling 
constants A k ^ are proportional to the probability den¬ 
sity of the electron at site k and the respective coupling 
strength A io with A x =Ay=—4f~—0.3iieV [140, 145 . 
The z-axis is chosen to be perpendicular to the graphene 
sheet. In contrast, the coupling strength of the (con¬ 
tact) hyperfine interaction in GaAs ranges from 74 /leV 
to 9 6fieV (depending on the specific isotope), thus ex¬ 
ceeding that of graphene by more than one order of mag¬ 
nitude [67 . The importance of the hyperfine interaction 
in graphene has been studied by comparing CVD-grown 
samples consisting entirely of 12 C and 13 C, respectively 
[146 . The experimental results indicate that in contrast 
to the cases of GaAs or carbon nanotubes, hyperfine in¬ 
teractions with 13 C nuclear spins have a negligible ef- 
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FIG. 8: We plot Eq. (18) for various parameters on a dou¬ 
bly logarithmic scale. The longitudinal hyperfine coupling 
strength in graphene is A z & 0.6 jueV [145 . The polarization of 
the nuclear spin system is assumed to be p— 0 for the solid lines 
and p— 0.9 for the dashed ones. For the thin black lines (left 
and bottom axes), we vary the relative 13 C abundance but fix 
the total number of atoms to iV=10 5 , which corresponds to 
a square QD with an approximate size of ^(50 nm) 2 . For the 
thick orange lines (right and top axes), we vary N and fix the 
relative 13 C abundance to its natural value, ^h^O.Ol. 


feet on spin transport relaxation times in CVD-grown 
graphene [67, 147, 148 . 

The interaction of an electron spin in a graphene QD 
with a certain number of nuclear spins, A13, will affect 
the time evolution of its transversal component (s + )(t), 
where s+=s x +is y . If Nis^>l, the nuclear spins form a 
bath and one can use the central limit theorem to obtain 

(s+)(t)=(s+}( 0)e T c [145 . In the presence of a perpen¬ 
dicular magnetic field B z ^> 2.6mT, the electron Zeeman 
energy is much greater than the transversal hyperfine 
coupling strength \A x ^ y \^0.3 fieV. Then, the character¬ 
istic decoherence time r c is given by 

2h N , 1CA 

where 0<p<l denotes the polarization of the nuclear spin 
system and N is the total number of carbon atoms, i.e., 
Ai3=0.01 for natural carbon [1451. Typical values for r c 
exceed 1 /is, see Fig. 8, In a further study, it was found 
that in the same situation, the longitudinal spin compo¬ 
nent ( s z )(t ) is conserved up to small corrections, which 
oscillate with a frequency determined by the hyperfine 
interaction [149 . If less than roughly 10 nuclear spins 
interact with the electron spin in the graphene QD, the 
nuclear spin system cannot be modeled as bath but al¬ 
lows for an exact treatment. Then, decoherence times 
depend on the exact position and orientation of the nu¬ 
clear spins, and typically he in the regime of 1ms [150 . 


VI. EXPERIMENTAL PROGRESS 

Spintronics with individual electron spins in semicon¬ 
ductors has evolved dramatically from basic concepts to 
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universal quantum control of spin states and quantum 
registers with up to four QDs [ 1 ,2,7,9, 38 , 39 , 58 , 63 , 64 . 
Isotopic purification enables coherence times beyond 1 s 
at low temperature in silicon and beyond 1 ms at room 
temperature in diamond [62 , 97 . 

Despite the challenge of efficient spin injection into 
graphene through tunneling barriers, remarkable im¬ 
provements in spin transport have been achieved over the 
past few years [12 , 69, 86 . In contrast to electron mo¬ 
bility, the spin lifetime seems to be rather insensitive to 
charged impurity scattering [70 and can exceed 1 ns at 
room temperature [71 , even for CVD-grown graphene 
[72 . The small intrinsic spin-orbit interaction makes 
it practically impossible to observe the spin Hall effect 
(SHE) in graphene [25 . But controlled hydrogenation 
can increase the spin-orbit interaction by three orders 
of magnitude and thus enable the observation of a SHE 
at room temperature and zero magnetic field [89 . Last 
year, a quantum spin Hall state has been observed in 
graphene with a large, tilted magnetic field [151 . 

Graphene QDs have been studied thoroughly 
w.r.t. their orbital spectrum [92—94, 99 , spin-filling 
sequence [152 , and charge relaxation times [153 j. More¬ 
over, high-quality bilayer QDs, encapsulated in hBN 
and equipped with electric gates, allow for electrostatic 
confinement as discussed in Subsec. Ill C. Small QDs 
in the 1 nm range can be created by electroburning and 
exhibit Coulomb blockade at room temperature due to 
addition energies as large as 1.6 eV [91 . Slightly larger 
graphene nanoflakes, epitaxially grown on Ir(lll) and 
intercalated with oxygen, are reported to exhibit a linear 
spectrum [154J. Despite these promising proceedings 
with graphene-based QDs, Pauli blockade — a basic 
step towards spintronics [155 and already achieved in 
carbon nanotubes [147, 148 — has (to our knowledge) 
not been observed in graphene, yet. 

Gapped armchair GNRs, Subsec. Ill B, and symmetri¬ 
cal graphene nanoflakes, Sec. IV, rely on precise edges. 
Atomically accurate armchair GNRs can be produced 
bottom-up, with widths up to 4nm, from aromatic pre¬ 
cursor molecules [156, 157 or top-down. One top-down 
technique that leads to precise edges is the unzipping 
of carbon nanotubes with an abruptly expanding nitro¬ 
gen gas inside the nanotubes [158 . Another possibil¬ 
ity is the electron-beam induced mechanical rupture of 
bulk graphene [159 . The latter method can be used 
to produce armchair edges as well as zigzag edges and 
can be applied in high vacuum, thus leading to mini¬ 
mally contaminated samples. Beyond this, clean struc¬ 
turing methods for graphene include nano-etching of sus¬ 
pended graphene [160 and the use of a silicon atom as a 
monatomic chisel [161 . 


VII. CONCLUSION AND OUTLOOK 

We have given an overview of spintronics with electron 
spin states in graphene quantum dots and nanoflakes. 
Long coherence times can be expected because of weak 
sources of decoherence, i.e., small coupling strengths of 
the spin-orbit and hyperfine interactions, and a marginal 
abundance of nuclear spins in carbon. We have explained 
spin relaxation via spin-orbit mediated coupling to acous¬ 
tic phonons and we have introduced a generic formula 
for the spin relaxation time T\ for graphene and other 
systems. Valley degeneracy and the lack of a bandgap 
in bulk graphene pose challenges to spin operations in 
graphene quantum dots. Specific systems where these 
challenges can be overcome have been reviewed and, if 
applicable, studies on their spin relaxation times have 
been discussed. Magnetism can be induced in graphene 
by several means. Here, we have focused on magnetism 
in graphene nanoflakes and in-situ tunability of antifer¬ 
romagnetism in symmetric nanoflakes with defects. We 
have summarized studies on the interaction of electron 
and nuclear spins and finally, we have pointed out recent 
experimental achievements relevant for spintronics with 
graphene quantum dots. 

We believe that graphene has the potential for a top- 
notch spintronics material that can outperform estab¬ 
lished materials like silicon or GaAs for specific applica¬ 
tions. The observation of the Pauli blockade represents 
a fundamental step towards spintronics with graphene 
quantum dots. In principle, the necessary architectures 
exist already but suffer from poor tunability of the source 
and drain barriers and inefficient spin injection. Simi¬ 
larly, the components required for electrostatic confine¬ 
ment of electrons in armchair GNR quantum dots are 
already available but have not been integrated into one 
device, yet. With view to the rapid progress of de¬ 
vice fabrication, we anticipate significant progress in 
these fields within the next 2-3 years. Ultimately, we 
see spintronics in graphene as a building block for 2D- 
heterostructured spintronics devices composed of (func¬ 
tionalized) graphene and other 2D materials like TMDCs 
and hBN. 
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